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Abstract-  Computed  Tomographic  Colonography  (CTC)  is  a 
minimally  invasive  method  that  allows  the  examination  of  the 
colon  wall  from  the  source  CT  sections  or  in  a  virtual 
environment.  The  primary  goal  of  CTC  is  the  detection  of 
colonic  polyps.  There  are  typically  two  data  sets  recorded  in 
supine  and  prone  positions.  An  important  step  in  polyp 
detection  is  the  anatomical  registration  of  these  datasets.  We 
developed  and  experimentally  validated  a  method  to  register 
complementary  supine  and  prone  CTC  datasets  anatomically, 
using  linear  stretching/shrinking  operations  on  the  medial  axis 
colonic  path  based  on  the  relative  path  geometries.  To  compare 
improvement  in  spatial  registration  of  supine  and  prone 
datasets,  a  radiologist  determined  5  unique  reference  points 
(RP)  in  the  registered  and  unregistered  datasets  from  each  of  5 
patients  by  viewing  supine  and  prone  data  simultaneously. 
Initial  results  suggest  that  our  algorithm  is  capable  of 
registering  supine  and  prone  CTC  data  anatomically  within  an 
approximate  range  of  ±13.2mm,  which  corresponds  to  0.8% 
error  with  respect  to  the  average  colon  length. 

Keywords  -  CTC,  Virtual  colonography,  registering  prone  and 
supine  data  for  CTC,  data  fusion,  anatomical  data  registration 

I.  Introduction 

Computed  Tomographic  Colonography  (CTC)  is  a  minimally 
invasive  method  for  the  examination  of  the  colon  using 
spiral/helical  CT  volume  data  [1],  The  standard  imaging 
process  consists  of  colon  cleansing  and  air-insufflation, 
followed  by  the  CT  imaging  of  the  abdomen  while  the  patient 
is  in  the  supine  position.  However,  insufficient  cleansing  and 
air-insufflation  often  cause  some  parts  of  the  colon  wall  to  be 
covered  with  water  [3-5].  Furthermore,  retained  solid  fecal 
material  can  mimic  a  polyp  in  CT  images  [2],  These 
limitations  are  largely  overcome  by  imaging  the  colon  a 
second  time  with  the  patient  prone.  The  radiologists  then 
compare  the  two  sets  of  images  simultaneously  to  assess  for 
movement  of  feces  and  to  examine  surfaces  hidden  in  supine 
position.  This  requires  careful  anatomic  alignment  between 
the  two  data  sets.  However,  as  the  patients  change  position, 
their  colons  deform  and  may  stretch  or  shrink  in  length. 
Currently,  anatomic  matching  the  two  data  sets  is  tedious  and 
time  consuming. 

In  this  study,  we  developed  and  evaluated  a  fully  automatic 
algorithm  for  the  registration  of  prone  and  supine  CTC  data. 
Such  an  algorithm  should  improve  the  radiologist’s  overall 
interpretation  efficiency  as  well  as  provide  a  basis  for 
combining  supine/prone  Computer-Aided  Detection  (CAD) 
results  automatically.  Our  approach  is  based  on  determining 
relatively  stationary  points  along  the  medial  axis  pathsof  the 
colon  (path)  for  both  positions,  then  matching  these  points  by 
simple  stretching  and/or  shrinking  of  either  the  supine  or 
prone  path.  The  algorithm  was  evaluated  on  a  set  of  5 
patients,  based  on  the  average  distance  along  the  path 
between  two  anatomically  identical  points  in  prone  and 
supine  data. 


Figure  1 .  X,  Y  and  Z  coordinates  along  the  supine  path  (Solid 
lines)  and  prone  path  (Dashed  lines).  Note  that  the  supine  path  is 
longer  than  the  prone  path. 


II.  Methodology 

A.  Background 

Colon  wall  segmentations  for  both  the  3D  prone  (pCT)  and 
supine  (sCT)  CT  data  were  done  using  a  fixed  intensity 
threshold  of  -700  HU.  The  paths  were  computed  using 
custom  software  and  sampled  at  1mm  intervals  [6].  The  X , 
y  and  z  coordinates  of  the  path  points  were  used  as  the 
inputs  to  the  registration  algorithm.  As  a  convention,  the  +  z 
axis  is  from  feet  to  head,  the  +  X  axis  is  from  the  left  side  of 
the  patient  to  the  right  side  and  the  +  y  axis  from  posterior 
to  anterior.  Distances  along  the  path  and  the  coordinates  are 
referenced  to  the  initial  path  point  at  the  anus. 

The  path  coordinates  of  sCT  (Xs,ys,Zs)  and  pCT 

( X  ,  y  , z p)  are  functions  of  the  distance  along  the  path, 
d  .  Fig.  1  shows  these  six  functions  for  a  single  patient.  An 
important  observation  at  this  point  is  that  x,y,Zs 
resemble  Xp,yp,Zp  ,  respectively  in  terms  of  their 

morphologies.  In  what  follows,  we  focus  on  the  Z  axis  for 
the  sake  of  simplicity,  though  the  algorithm  can  be  similarly 
applied  to  any  individual  (or  combination  of)  axis  (axes). 
This  latter  approach  is  addressed  later. 

The  morphological  similarity  between  zs  and  Z. p  can  be 

seen  as  a  similarity  between  the  number  of  local  extrema 
(LoP)  and  their  order  in  terms  of  their  z  values.  In  other 
words,  if  we  sort  the  LoP’s  with  respect  to  their  z  values,  the 
LoP’s  of  Zs  and  Z,  p  with  equal  or  close  indices  are  more 

likely  to  correspond  to  the  same  anatomical  points  along  the 
paths.  In  anatomic  terms,  LoP’s  correspond  to  the  deflection 
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Figure  2.  The  colon  wall  extracted  from  3D  CT  data  by  simple 
thresholding  and  the  central  axis  path  extending  from  the  anus  to 
the  cecum. 

points  along  the  path.  Despite  of  the  fact  that  these  deflection 
points  move  within  the  3D  volume  (the  patients’  abdomen), 
they  remain  to  be  deflection  points  with  small  shifts  in  their 
locations  with  respect  to  d  .  Furthermore,  LoP’s  with  high  z 
values  are  the  points  furthest  away  from  the  anus. 
Anatomically,  these  points  typically  correspond  with  the 
hepatic  and  splenic  flexures,  locations  that  are  relatively  fixed 
in  location.  These  flexures  mark  the  beginning  and  ending  of 
the  much  more  mobile  transverse  colon  (see  Fig.  2).  (LoP’s 
along  the  other  axes  correspond  with  different  anatomical 

points.)  Our  method  is  based  on  coupling  LoP’s  of  zs  and 
Zp  and  matching  their  positions  with  respect  to  d  by  simple 
segmental  stretching  and/or  shrinking  of  one  of  the  two  paths. 

B.  Algorithm 

The  first  step  is  baseline  adjustment.  We  assume  that  the 
beginning  and  the  ending  points  of  the  two  paths  (spPt  and 
prPt)  correspond  to  the  same  anatomical  points,  the  anus  and 
the  tip  of  the  cecum.  Let  the  path  length  functions  of  spPt 
and  prPt  be  Ds  (z  )  =  i,  i  €E  [1,  Ns  ]  and 

D  (l)  =  i,  i  €E  [1,  N  ] ,  respectively.  A  linear 
transformation,  T  ,  is  applied  such  that 

Dp(l)  =  Ds(l)  =  T(Ds(l))  (1) 

Dp(Np)  =  Ds(Ns)  =  T(Ds(Ns ))  (2) 

This  is  equivalent  to  taking  the  first  and  the  last  points  of  the 
two  paths  as  two  pairs  of  coupled  LoP’s  of  spPt  and  prPt  and 
matching  them  with  respect  to  d  .  The  z  values  along  the 
spPt  are  also  changed  linearly  such  that  Zs  (1)  =  Zp  (1)  and 

Zs  (Ns )  =  Zp  (./V  )  .  Note  that  the  linear  transformation, 
T  ,  is  the  mapping  between  spPt  and  prPt.  It  will  be  modified 


Figure  3.  Z  coordinates  versus  path  length  of  supine  (  Z  s )  and  prone 

(Z p)  data  after  baseline  adjustment.  The  *  indicate  the  LoP_s’s  and 
LoP_p’s 

in  the  subsequent  steps  in  a  nonlinear  fashion  to  yield  to  final 
piecewise  linear  mapping  function. 

The  second  step  is  to  determine  the  LoP’s.  The  local 

extrema  of  ZJ(Z)S)  and  zp{Dp)  are  determined  using 

numerical  derivation  and  neighboring  extrema  are  clustered 
such  that  the  minimum  distance  along  the  path  between  two 
consecutive  clusters  is  15mm.  The  mean  location  along  the 
path  and  the  mean  z  value  for  spPt  and  prPt  (LoP_s  and 
LoP_p  respectively)  are  used.  LoP_s’s  and  LoP_p’s,  thus 
found,  are  ordered  with  respect  to  their  z  values  and 
indexed.  Fig.  3  shows  the  LoP_s’s  and  LoP_p’s  for  Z  and 


The  next  step  is  to  find  the  LoP_s,-LoP_p  pair  most  likely 
to  correspond  to  the  same  anatomical  point.  The  decision 
criteria  is  based  on  the  differences  between  their  z  values 
and  their  locations  along  the  path: 

If 

i,  j  e  {1,2}  and  (3) 


D  ( LoP_s(i))-D  {LoP _  p{  j)\  <  100mm  (4) 

s  p 


K  <  D  ,D  <  L 
s  p 


(5) 


where  i  and  j  are  indices  assigned  to  LoP’s  in  the  region  of 
interest  (ROI,  Region  Of  Interest,  is  defined  by  (5))  only, 
with  respect  to  their  z  values.  Condition  (3)  forces  the 
algorithm  to  match  LoP’s  starting  from  the  ones  furthest 
away  from  the  anus,  i.e.,  the  ones  most  likely  to  correspond 
to  anatomical  landmarks  (like  the  hepatic  and  splenic 
flexures).  Condition  (4)  is  a  constraint  on  the  maximum 
stretching/shrinking  allowed.  Condition  (5)  defines  the 
boundaries  of  ROI  on  which  the  registration  is  applied. 

Initially,  K  =  Dp  (N  p  )  =  Ds  (N. )  and 
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1. Do  baseline  correction 

2.  Find  LoP’s 

3. Initialize  recursion_level=0  and  K  and  L  to  set  ROI 
4. Sort  and  index  LoP’s  in  ROI 
5  .If  recursion_level<max_recursion_level 
5.1. Match  a  pair  ofLoP_s  and  LoP _p  based  on  (3), (4), (5) 
5. 2. If  no  match,  go  to  step  6 

5. 3.  Update  D .  by  applying  T{  ) ,  as  in  (6) 

5. 4. Increase  recursion_level  by  1 
5.5.  K0  =  K,  L0  =  L 

5. 6. Set  K  =  Ds  ( LoP  _  s(i ))  =  Dp  ( LoP  _  p(j )) 

5. 7.  Go  to  step  4  ( Process  right-hand  side) 

5. 8.  Set 

K  =  K0,  L  =  Ds  ( LoP  _  s(i))  =  Dp  ( LoP  _  p(j )) 

5. 9.  Go  to  step  4  (Process  left-hand  side) 

6. Output  Z)s. 


Figure  4.  The  recursive  registration  algorithm 

L  =  D p  (1)  =  Ds  (1)  .  Then  a  linear  stretching/shrinking 
operation  is  applied  to  Ds  such  that 

T(D,  ( LoP  _  s(i)))  =  Dp  ( LoP  _  p(j ))  (6) 

The  algorithm  matches  the  right-  and  left-  hand  sides  of  the 
matched  point  recursively,  for  a  pre-specified  level  of 
recursions  or  until  it  could  not  find  any  pair  of  LoP’s  that 
satisfy  the  above  criteria.  Fig.  4  summarizes  the  algorithm. 

C.  Iteration  and  Parameters 

The  algorithm  described  above  for  z  axis,  can  be  similarly 
applied  to  X  and  y  axes  exactly  in  the  same  way.  The  order 
in  which  the  axes  are  considered  is  a  parameter  of  the 
algorithm.  The  entire  process  can  be  iterated  until  it 

converges  to  a  final  D  or  for  a  pre-specified  number  of 
iterations. 

Table  1  provides  a  list  of  the  parameters  of  the  algorithm. 
The  analysis  in  the  following  section  will  be  based  on 
variations  of  these  parameters. 

D.  Initial  Clinical  Evaluation 

We  evaluated  the  improvement  in  spatial  registration  of 
identical  anatomy  in  supine  and  prone  data  sets  from  5 
patients  (all  male,  age  64-77,  mean  age  70).  CT  data  were 
acquired  in  each  position  after  colon  cleansing  and  air- 
insufflation,  using  a  single-detector  CT  system  (GE  Medical 
Systems,  Milwaukee,  WI).  Imaging  parameters  were  3mm 
collimation,  pitch  1. 5-2.0,  1.5mm  reconstruction  interval. 


Table i 

THE  PARAMETERS  OF  THE  ALGORITH 


Parameter 

Definition 

Maxlter 

Maximum  number  of  iterations 

Axes 

Axes  to  be  used  and  their  order 

MaxRec 

Maximum  recursion  level 

ClstDLoP 

Minimum  d  between  neighboring 

groups  of  local  extrema 

MaxDLoP 

Maximum  d  between  matched  LoP_s 
and  LoP_p.  Different  for  X,  y,  Z  ■ 

120  KVp,  200  MAs.  A  radiologist  determined  5  unique 
reference  points  (RP)  in  each  of  5  patients  by  viewing  supine 
and  prone  data  simultaneously.  Reference  points  included 
polyps,  diverticulae,  unique  haustral  folds,  or  the  ileocecal 
valve.  We  computed  the  distance  in  terms  of  d  between 
anatomically  identical  RPs  in  supine  and  prone  data  both 
before  and  after  path  registration.  The  change  in  the  mean 
distance  between  pre-  and  post-registration  RP  positions  was 
computed  as  a  percentage  of  the  pre -registration  distance.  In 
the  following,  we  report  the  mean  of  these  percentage  values 
( Ad )  over  all  patients  in  the  dataset.  Ad  provides  an 
overall  understanding  of  the  performance  of  the  algorithm 
and  is  used  to  assess  different  parameter  values.  The 
parameter  ranges  used  are:  Maxlter=[l  2],  Axes=[{  z  }, 
{  Z,x,y },  {z,x},  {z,y}],  MaxRec=[l  2  3], 

ClstDLoP=[10  15  ...  40],  MaxDLoP=[(42,42,70),  (48,48,80), 
...,  (72,72,120)].  Each  parameter  was  evaluated  separately, 
keeping  the  others  constant. 

III.  Results 

To  assess  each  parameter,  we  assigned  a  value  to  that 
parameter  from  the  sets  given  above,  held  it  constant  and 
computed  Ad  for  all  combinations  of  the  remaining 
parameters.  This  resulted  in  504,  252,  336,  144,  and  168 
Ad ’s  for  a  single  value  of  Maxlter,  Axes,  MaxRec, 
ClstDLoP,  and  MaxDLoP,  respectively.  The  respective  mean 
values  of  Ad ’s  were  [-52.3%,  -52.5%],  [-52.4%,  -52.0%,  - 
53.3%,  -52.0%],  [-52.5%,  -52.4%,  -52.4%],  [-35.6%,  -56.4%, 
-53.8%,  -51.2%,  -55.3%,  -57.2%,  -57.5%],  [-44.4%,  -47.5%, 
-47.3%,  -59.5%,  -60.4%,  -55.5%].  Note  that  more  negative 
the  mean  Ad  ,  the  better  registered  the  RPs.  We  assessed  for 
significant  differences  between  Ad  ’s  between  values  of  a 
single  parameter  with  the  Student’s  t-test.  All  p-values  were 
less  than  0.05  except  in  the  following  cases:  (i)  between 
different  MaxRec  values  and,  (ii)  between  using  axes  {  z  }, 
{  z,X,y  }  and  {  Z,y  }■ 

Based  on  these  results,  we  decided  to  use  the  following  set 
of  parameters  for  path  registration  in  this  data  set,  with  our 
rationale  explained  following  each  parameter:  i)  Maxlter=l: 
The  difference  in  mean  Ad ’s  was  small,  ii)  Axes={  z,X  }: 
There  is  a  significant  improvement  of  approximately  1%  with 
respect  to  the  other  choices,  iii)  MaxRec=l:  Higher  values 
did  not  improve  the  performance  significantly,  iv) 
ClstDLoP=15:  This  corresponded  with  a  local  minimum 
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Figure  5.  Pre-  and  post-registration  distances  along  the  path  between 
25  pairs  of  RPs  in  supine  and  prone  CTC  data. 


significantly  different  that  the  others.  The  lower  mean  A  d 
values  for  ClstDLoP=35  and  40  were  not  chosen  as  they  did 
not  represent  a  local  optimum  point  and  because  large 
ClstDLoP  values  result  in  large  clusters  that  may  localize 
LoP’s  with  sufficient  resolution,  v)  MaxDLoP=[66  66  100]: 
This  was  a  local  minimum  significantly  different  from  the 
other  values. 

The  mean  (±s.d.)  Ad  with  the  optimal  parameter  settings 
was  (-67.0±25.8)%.  The  mean  (±s.d.)  pre -registration 
distance  of  51.0142. 8mm  was  decreased  to  14.2+1 1.3mm. 
The  maximum  pre -registration  distance  of  151.0mm  was 
decreased  to  22.5mm.  Fig. 5  summarizes  the  results.  The 
point-to-point  distance  was  increased  by  registration  in  only  4 
out  of  25  RPs  (5  RP  per  patient)  (from  14.3+1 1.9mm  to 
23. 0+20. 5mm).  These  are  the  points  above  the  y=x  line  in 
Fig.  5. 

IV.  Discussion 

Fig.  5  shows  every  RP  in  the  dataset  with  the  line  fitted  in 
the  least-squares  sense.  It  suggests  that  our  method  can 
register  any  two  points  along  the  colon  path  with  an 
approximate  average  error  of  14.2mm,  irrespective  of  the  pre¬ 
registration  distance  between  these  points. 

These  initial  results  support  the  validity  of  our  assumption 
that  the  colon  undergoes  apparent  linear  stretching  and 
shrinking  with  changes  in  patient  position.  There  may  also 
occur  twisting  of  the  colon  around  the  colon  path,  though  we 
do  not  yet  evaluate  for  this  effect.  Such  evaluation  will 
require  a  more  detailed  point-wise  registration  of  the  colon 
walls  in  different  positions. 

One  of  our  results  presented  above  show  that  the  larger  the 
variation  in  the  coordinate  values,  the  more  the  information 
content.  In  support  of  this,  we  found  that  the  Z  and  X  axes 
are  the  most  important  axes.  This  makes  sense  as  the  colon 
more  or  less  lies  on  the  ( z ,  x)  plane  (coronal  plane). 


V.  Conclusion 

We  have  developed  and  experimentally  evaluated  an 
algorithm  to  register  supine  and  prone  CTC  data  with  respect 
to  positions  along  the  central  path  of  the  colon.  Possible 
second  order  effects  such  as  twisting  of  the  colon  wall  around 
the  path  is  the  subject  of  future  research.  These  initial  results 
suggest  that  our  algorithm  is  capable  of  registering  supine 
and  prone  CTC  data  anatomically  within  an  approximate 
average  error  of  14.2mm,  which  corresponds  to  0.9%  error 
with  respect  to  the  average  colon  length.  This  is  a  significant 
improvement  over  the  pre -registration  distance  of  51mm. 
From  a  clinical  perspective,  this  suggests  that  a  given 
anatomic  point  in  a  complementary  supine  or  prone  dataset 
can  be  identified  with  a  ~3  cm  zone  compared  with  a  -10  cm 
range  using  non-registered  data.  With  further  validation  and 
refinement,  our  path-registration  algorithm  will  allow  for 
more  efficient  simultaneous  data  viewing  by  radiologists  and 
enable  accurate  anatomic  correlation  of  computer-aided 
detection  results  between  patient  positions. 
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